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1 Introduction 



The connection between astrophysics and laboratory experiments is often 
made in the context of atomic physics, and in particular spectroscopic stud- 
ies. These generally concentrate on the micro-physics of astrophysical plas- 
mas, providing important data on cross-sections, opacities, grain chemistry, 
and other processes that are vital to the modelling and interpretation of as- 
tronomical observations. Indeed, detailed laboratory opacity measurements of 
Q iron absorption lines have greatly improved models of pulsation periods in 

^ Cepheid variables [TTj. 

Recreating large-scale astrophysical phenomena in the laboratory has also 
inspired scientists. A pioneer in the use of laboratory experiments to study 
space physics was the Norwegian scientist Kristian Birkeland (1867-1917), who 
used a terrella[^ and gas discharges to investigate auroral phenomena (22j . In- 
>^ deed, as early as the 17th century a terrella was used by William Gilbert to 

study the Earth's magnetism f3D] . More recently, the idea that physical prob- 
lems in the cosmos may be elucidated by more mundane earth-based events, 
such as fluid flows and nowadays laboratory plasmas, was taken up in a series 
of meetings held just after World War II, with the aim: " To bring together 
workers from astrophysics and from aerodynamics; ...to consider which devel- 
^\ opments in fluid mechanics may be applicable to astrophysical problems, and to 

arrive at a formulation of these problems in such a way that mathematicians 
^ and fluid mechanics people may find a way of attack" ^55,. Topics of discus- 

sion included accretion, gas and dust in the interstellar medium, shocks, and 
turbulence; the latter being very popular with those working in fluid mechan- 
ics. An indication that this was taken rather seriously can be gleaned from 
the list of scientist who took part in the meetings. These included the likes of 
Heisenberg, Sciama, Bondi, Hoyle, Bok, Burgers, Seaton and Cowling. At the 
meeting held in Cambridge in 1954, Chandrasekar and Fermi contributed a 



^ A terrella is a small (tens of cm) magnetized model Earth. 
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preliminary paper but did not actually attend; while Hubble, although listed 
as participant, had died in fact nine months earlieij^ 

In the last ten years or so, experiments on pulsed-power facilities (z- 
pinches) and high-power lasers are leading the way in studying astrophysical 
phenomena in the laboratory. Emerging areas of research have been aimed 
at producing complex dynamical phenomena, such as compressible hydrody- 
namic mixing, hypersonic jets, shock physics, radiation hydrodynamics and 
photoionised plasmas, to name a few. These can help to understand the physics 
of phenomena associated with wide range of astrophysical objects, including 
protostellar and AGNs jets, supernovae explosions and the subsequent genera- 
tion of remnants, and photoevaporated molecular clouds. Here we will restrict 
the discussion to the jets produced on z-pinch facilities and in particular to 
the work performed on the MAGPIE generator at Imperial College. For a 
broader and more detailed discussion of laboratory astrophysics on lasers and 
z-pinches we refer the reader to the review by Remington et al. 2006f46| and 
the book by Drake 2006 



2 Plasma conditions in z-pinch and laser experiments 

Typical plasmas produced on z-pinch and laser facilities have pressures of ^ 
Mbar, corresponding to energy densities ~ 10^^ erg cm~^, at a fraction of 
solid density. An overview of the plasma conditions attainable on experimen- 
tal installations, together with some of those found in space is given in Fig. [l] 
Z-pinch facilities rely on stored electrical energy (hundreds of kJ) to deliver 
large currents (~ of a few Mega Amperes) over a short time (~ 100 — 1000 
ns) to a "load" usually consisting of a gas or thin metallic wires. These fa- 
cilities typically produce volumes of plasma ~ 1 cm~^ (for a review see |48j'). 
Laser facilities instead rely on focusing onto a solid or gaseous target, single or 
multiple high-power laser beams. These concentrate several kJ of energy, over 
timescales ~ pico- to nano-second, into plasma volumes of ~ 1 mm"'^. The 
future arrival of the Laser Megajoule (LMJ) facility in France and the Na- 
tional Ignition Facility (NIF) in the USA, will produce fusion plasmas under 
conditions similar to stellar interiors. 

For the present discussion (i.e. jets) it is more interesting to look at the dy- 
namical conditions that can be obtained in the laboratory. When the energy 
available on z-pinches and lasers is partially converted into kinetic energy, 
it can generate hypersonic (Mach numbers Af > 5), radiatively cooled flows 
with characteristic velocities of the order of 100 - 1000 km s~^. These flows 
can include dynamically important magnetic fields, ~ several 10^ Gauss, and 
have a large range of plasma-/? (the ratio of thermal to magnetic pressure), 
1 >> (3 >> 1. In such cases, our inability to obtain the adequate astrophysical 



As a final piece of trivia, Hubble's funeral was never held and the location of his 
resting body never disclosed. 
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plasma conditions (see Fig. [T| may be overcome by producing scaled "condi- 
tions" of the phenomena of interest. These are discussed in the next section. 



3 Relating laboratory and astrophysical phenomena 

The framework to relate experiments to astrophysical phenomena, within 
magnetohydrodynamics, was developed in a series of papers by Ryutov et 
al. 1371 Hni Uni [51] . A general discussion of hydrodynamic scaling can also be 
found in [SB], and in the present volume it is reviewed in detail by Cocker. An 
interesting duality between imploding and exploding systems is given in [43] . 
For coUisionless plasmas the scaling conditions are descibed in [20 . 

Here we will qualitatively introduce some of the ideas behind scaling and 
shall be concerned with systems which are to a good degree magnetofluids. 
As a specific example we take the MHD jet launching, which will be later 
discussed in more details in the context of laboratory experiments. Fig. [2] 
shows schematically the comparison of the "environment" producing astro- 
physical and laboratory jets. The details of the experimental set-up depicted 
are not important at this stage, it is sufficient to say that the physical pro- 
cesses leading to the distributions of magnetic fields, plasma density, pressure 
and velocity inside the "modelling box" (MB) are completely different for the 
two systems. However because we are interested in studying experimentally 
the collimation and launching of the jet in the MB, we are concerned only with 
the differences that may be present there. The physical model to be employed 
in the MB is that of ideal MHD (for an introduction to MHD see gl])- While 
this is a good approximation for the astrophysical case, it is not usually the 
case for the experiments, which have to deal with non-ideal effects arising for 
example from thermal conduction, viscosity and resistivity. Nevertheless the 
experiments can be designed so that the ideal MHD approximation is valid at 
least inside the MB, thus making the set of equations describing the plasma 
evolution in the two MB the same. For scaling purposes this is however not 
sufficient, it is also fundamental for the initial (and boundary) conditions in 
the two MB to be geometrically similar. Meaning that at some arbitrary time, 
the initial spatial distributions of all physical quantities (density p, velocity 
V, pressure P and magnetic field B) are the same. Finally, the actual scal- 
ing factors relating the physical variable in the two systems and their initial 
distributions should obey the following constraints: 



where the tilde denotes characteristic (dimensional) values for the laboratory 
(L) and astrophysical (A) plasmas. If the aforementioned conditions are met, 
namely the ideal MHD equations hold, the initial conditions are geometrically 
similar and the constraints on the scaling factors hold, then the evolution in 
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the modelling boxes, from some initial time and over the scaled time span, 
will be indistinguishable. The addition of extra physical processes, for example 
radiation, increases the constraints to be satisfied and hence the difficulty in 
obtaining an exact scaling (see for example 

Two points need to be clarified. The first regards obtaining geometrically 
similar initial conditions in the laboratory. This is clearly very difficult and it 
has been done exactly in only a few cases [45]. The question also concerns which 
initial condition should be used. Obviously observations cannot be expected to 
provide three-dimensional distributions of all the quantities of interest , if they 
did the problem would be almost solved! Then, with observational constraints 
in mind, we resort to models, which are in general as difficult to reproduce. 
The usefulness of laboratory astrophysics however goes beyond strict scaling. 
By using "reasonable" initial conditions, as it is usually done in numerical 
simulations, we can study in detail complex astrophysical phenomena in a re- 
peatable and accessible manner. In addition the possibility of validating astro- 
physical codes using laboratory data is also very important. The second point 
is on the accuracy of the approximation of ideal MHD for laboratory plasmas, 
and more specifically the ones produced on z-pinch installations. Neglecting 
the effects associated with viscosity, thermal conduction and resistivity (the 
ideal MHD approximation) means that the Reynolds (Re), Peclet (Pe) and 
magnetic Reynolds (Re^) numbers are Re, Pe, Re^ >> 1- While for YSO 
jets these dimensionless numbers are very large (see Table [T]), for laboratory 
jets they are generally many orders of magnitude smaller, thus calling into 
question the validity of the ideal MHD approximation. However the values 
obtained in the experiments should be compared to astrophysical numerical 
simulations, where because of the finite accuracy of the numerical schemes 
employed, the ideal MHD approximation also breaks down. Unless included 
explicitly in the simulations, viscous and resistive dissipation, and thermal 
conduction occur at the grid level through unphysical numerical truncation 
errors. Existing ideal MHD simulations of jets have dimensionless numbers in 
the range 50 — 1000, well within the reach of laboratory experiments. In some 
case, as for the Reynolds number, the laboratory values can be orders of mag- 
nitude larger. Thus laboratory astrophysics experiments can be very useful in 
studying complex and intrinsically non-linear three-dimensional phenomena, 
in regimes where accurate numerical simulations are particularly difficult to 
perform, and thus it represents a powerful complementary tool to astrophysi- 
cal modelling. As a final justification: "Real experiments are also irreplaceable 
in providing new insights into subtle physics issues and in stirring the creative 
imagination of scientists" [50] . 

4 Young stellar jets from Z-pinch machines 

In this section we shall discuss the z-pinch studies of astrophysical jets. It is 
worth mentioning that the time for an experiment to reach "maturity" : from 
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the conception, design and first experiments to the modeUing, data analysis 
and then the apphcation to astrophysical models is of the order of 3 to 5 years, 
depending on its complexity. Considering that most of the experiments, both 
on laser and on z-pinches, are less than 10 years old, it should not come as a 
surprise that in many instances the applications to astrophysics are not yet 
fully developed or very clear. Although we shall restrict our discussion to jets, 
there are other areas of laboratory astrophysics research on z-pinches that 
are well developed, such as equations of state studies for planetary interiors 
and photoionised plasmas relevant to accretion disks around compact objects. 
These are thoroughly reviewed in |46] and we refer the interested reader to 
that paper. 

4.1 Hydrodynamic jets from conical wire arrays 

Conically converging flows were investigated in astrophysics as a possible 
mechanism for converting wide angle winds into coUimated jets. Such models 
do not require magnetic fields, at least to coUimate and launch the jet, and 
rely instead on purely hydrodynamic means [52l [6]. Within this framework, 
and with the aim of producing hydrodynamic jets to be used for interaction 
studies, a series of experiments [3TJ [TD] were developed on the z-pinch gener- 
ator MAGPIE. The schematic of the experimental configuration is shown in 
Fig. [Sj It consists of a conical array of micron-sized metallic wires driven by 
a current of 1 MA rising to its peak value in 240 ns. The basic mechanism 
of plasma formation in wire arrays is the following: resistive heating rapidly 
converts the wires into a heterogeneous structure consisting of a cold (< 1 
eV) and dense liquid- vapour core, surrounded by a relatively hot (10 — 20 eV) 
and low density (~10^^ cm"-^) plasma. Most of the current fiows in the latter, 
where the resistivity is lower, which undergoes acceleration by the J x B force 
toward the array axis. These streams of plasma have characteristic velocities 
of ~ 100 — 150 km and corresponding Mach numbers M ^ 5. The wire 
cores act as a reservoir of plasma, replenishing the streams during the entire 
duration of the experiment (several hundred ns). The converging plasma is 
virtually magnetically field-free and the interaction on axis is hydrodynamic 
in character. The collision produces a standing conical shock where part of the 
kinetic energy of the streams is thermalized. However it is important to note 
that the plasma streams are not perpendicular to the surface of the conical 
shock. Thus the component of the velocity parallel to the shock is continuous 
across the shock and the flow is redirected upwards into a jet. Typical jet 
velocities are ~ 100 — 200 km s~^ and hypersonic jets with Af > 10 can be 
produced by this mechanism. The jet coUimation and Mach numbers depend 
predominantly on the level of radiation cooling in the plasma, which can be 
changed experimentally by varying the wire material (Al, Fe, W and so on). In- 
creasing the atomic number of the wires increases the rate of cooling from the 
plasma, lowers its temperature and leads to the formation of more collimated 
jets (with higher Mach numbers) [lOl |3TJ |32] . These jets are used to study the 
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propagation and interaction with an ambient medium, which are described in 
sec. |4.3[ The characteristic conditions and dimensionless parameters obtained 
are shown in Table [TJ 

It is possible to design and modify the experiments to include additional 
physics, such as dynamically dominant magnetic fields and rotation. Indeed, 
for accretion onto the forming star to occur, angular momentum needs to 
be removed from the in-falling material. In combination with the processes 
present in the accretion disk, such as MHD instabilities and turbulent trans- 
port, jets and winds can also remove a considerable fraction of the excess 
angular momentum from the accreting flow [16, . One of the obvious implica- 
tions is that jet will be rotating and some confirmation has arrived with recent 
observations of rotation in a number of YSO jets [H]. Supersonically rotating 
laboratory jets and fiows of astrophysical relevance were recently produced 
for the first time |5] using a variant of the conical wire array. Rotation in the 
flow is accomplished by slightly twisting the wires in the azimuthal direction. 
This results in a poloidal magnetic fleld and an azimuthal component of the 
Lorentz force, giving a non-zero torque on the plasma streams (Fig. |4]). The 
level of angular momentum introduced in the system can be controlled by 
changing the twist angle and in general the jets ejected have rotation veloci- 
ties ~ 100 — 200 km s^^, corresponding to ~ 0.1 — 0.2 of the jet propagation 
velocity. One of the applications of these proof-of-principle experiments will 
be to study the effects of rotation on the propagation of jets and on the growth 
of the Ray leigh- Taylor instability in curved jets (see section 4.3 1 



4.2 Magnetohydrodynamic jets 

Protostellar (and galactic jets) are thought to be powered by the combination 
of rotation and magnetic fields, which extract the rotational energy from an 
accreting system and create magnetic stresses which accelerate and coUimate 
the fiow (see the lectures notes |23l[53] . Depending on the details of the models, 
the winding of an initially poloidal magnetic field results in a fiow pattern 
dominated by a toroidal field. A similar situation is also attained when the 
foot-points of a field line, connecting the disc to a central compact object or 
connecting different parts of a disc, rotate with different angular velocities. 
In such cases, the relative angular displacement of the foot-points causes one 
of them to move ahead of the other, and the field loop to twist. The induced 
toroidal component results in an increase of the magnetic pressure which 
drives the expansion of the loop itself |Mj . In the magnetic tower scenario [571 
I36j . the outcome is a magnetic cavity consisting of a highly wound up toroidal 
field which accelerates the fiow. In this case, the presence of an external plasma 
medium was shown to be necessary to confine the magnetic cavity, which 
would otherwise splay out to infinity within a few rotations [35j. The basic 
picture of magnetic tower evolution has also been confirmed numerically by 
several authors [251 Hi SD il ISH] ■ 
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Magnetically driven jets from radial wire arrays 

The study of magnetically coUimated and accelerated jets on z-pinches was 
developed in the last few years [HI I3S1 El [H] using a modified wire array 
configuration. The basic astrophysical mechanism studied in the experiments 
is the interaction of a toroidal magnetic field with a plasma ambient medium, 
leading to the formation of jets and magnetic "bubbles" . The schematic of 
the experimental set-up, a radial wire array, is shown in Fig. [5] The forma- 
tion of plasma is similar to that discussed in conical wire arrays, however 
the plasma is now accelerated vertically filling the space (few cm) above the 
array. Below the wires there is only a toroidal magnetic field. The formation 
of the jet and its time evolution is shown in Fig. [6] The initial formation of 
the magnetic cavity and jet occurs at the time when the magnetic pressure 
is large enough to break through the wires. This occurs only over a small 
region close to the central electrode, where the toroidal magnetic field Bq is 
strongest. The results show the system evolving into a structure consisting 
of an approximately cylindrical magnetic cavity with an embedded jet on its 
axis confined by the magnetic "pinching" force. A shell of swept-up plasma 
surrounds and partially confines the magnetic bubble. The subsequent evolu- 
tion is dominated by current-driven instabilities and the development of the 
asymmetric "kink" mode (m — 1) which leads to a distortion of the jet and a 
re-arrangement of the magnetic field. In Fig. [7]i, the magnetic field lines can 
be seen to twist inside the jet, an effect caused by the instability which turns 
toroidal into poloidal magnetic flux. The end result of the instabilities how- 
ever is not to destroy the jet, but to produce an inhomogeneous or "knotty" 
jet, shown in Fig.[7]3-c. The resulting jet has typical super-fast-magnetosonic 
Mach numbers in excess of 5, it is kinetically dominated and its opening angle 
< 20°. 

The relatively simple initial conditions implemented experimentally pro- 
duce a very complex and rich dynamics which share many important features 
with astrophysical models. One important example is the presence of a en- 
velope surrounding the magnetic cavity and confining it. Although this is 
discussed in the astrophysical literature |37l|Sl], it has so far only been ob- 
served in a laboratory experiment. The stability and dynamics of the envelope, 
which determine the coUimation of the cavity itself, can thus be directly stud- 
ied in the laboratory before astronomical observations may become available. 
Finally, it is worth pointing out that while two-dimensional, axisymmetric 
MHD simulations reproduce very well the experimental results, up to the de- 
velopment of the non-asymmetric current-driven instabilities. There are fun- 
damental differences in the long-term evolution of the system, which can only 
be reproduced by fully three-dimensional simulations. 

Episodic ejection of magnetic bubbles and jets 

Protostellar jets are characterized by the presence of knots and multiple bow- 
shocks [2F, tracing their propagation. These are often interpreted as internal 
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shocks driven by relatively small perturbations in a steady ejection process, 
and which occur on typical time-scale between ~ 5 — 20 years. For example 
in [TH] it was shown that the temporal variability of the jet velocity may 
be associated with a time-varying stellar magnetic field. Episodic jet ejection 
behaviour may also be associated with variation in the accretion rates or an 
inflating stellar magnetosphere. 

Recent experiments have studied for the first time the episodic ejection of 
magnetic bubbles and jets, and its effects on the overall propagation of the 
outfiow[T5]. The experimental set-up is similar to that shown in Fig. |4] how- 
ever the wires are replaced by a 6 /xm thick metallic foil (usually aluminium) . 
A 3D MHD simulation of the experiments is shown in Fig. [8] The evolut ion 
of the first bubble is similar to that of radial wire arrays. However the total 
mass in the plasma source, as a function of radius, is larger for a foil than for 
a radial wire array. Thus after the first magnetic cavity and jet are formed, 
there is a larger quantity of plasma available to refill the "gap" between the 
central electrode and the left over foil; for example the presence of this gap 
is visible in|5]for radial arrays, and it is produced by the magnetic field pres- 
sure breaking through the wires or foil. Once the gap is refilled with plasma, 
the currents can fiow once again across the base of the magnetic cavity, thus 
re-establishing the initial configuration. When the magnetic pressure is large 
enough to break through this newly deposited mass, a new jet/bubble ejec- 
tion cycle can begin. Typical fiow velocities observed are ~ 100 — 400 km 
s~^, the simulated sonic and the alfvenic Mach numbers in the jet, defined 
as the ratios of the fiow speed to the sound and Alfven speed respectively, 
are Ms ~ Ma 3 — 10. The resulting fiow is heterogeneous and clumpy, and 
it is injected into a long lasting and well coUimated channel made of nested 
cavities. Each jet/outfiow episode propagates, interacts and substantially al- 
ters the surrounding environment by injecting mass, momentum, energy and 
magnetic flux into it. An important aspect of the episodic ejection process 
is, broadly speaking, its self-coUimation. Since the initial ambient medium is 
swept away after a few ejections, newly formed magnetic cavities are confined 
solely by the environment left by earlier episodes, thus making the coUimation 
process insensitive to the initial ambient conditions. An experimental image 
of the evolution of the system is shown in Figj9] In the magnetic cavities 
Rem > 100, and each bubble expands with its own "frozen- in" magnetic fiux; 
in the experiments this is confirmed by the magnetic probe measurements of 
the trapped magnetic field at the outer edge of the bubbles, B ^ 1 — b kG. The 
coUimation is then determined not only by the pressure of the left-over plasma 
but also by the pressure of the tangled magnetic field trapped in the bubbles, 
where the plasma-/? is in the range 0.1 < /3 < 1. A high level of symmetry 
is maintained after many ejections (5 in the current experiments), the num- 
ber being limited only by the duration of the current pulse delivered by the 
generator. Overall, the experiments demonstrate that magnetic acceleration 
and coUimation, occurring within a framework of strongly episodic outflow 
activity, can be effective in producing well coUimated and heterogeneous jets. 
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By drawing a parallel with the dynamics observed in the experiments, one 
can gain useful insights and a qualitative view of the possible evolution of as- 
trophysical jets. In the experiments there are two time-scales which determine 
the magnetic bubbles/jets development: the current-driven (CD) instability 
time-scale t/ and the episodic bubble ejection time-scale tb- For conditions 
applicable to the formation region of protostellar jets |26] we can estimate 
the growth time of the CD kink mode as the Alfven crossing time r/ ~ 1 
year; corresponding to a few nanoseconds in the experiments. The second 
time-scale is linked to the temporal variability of the Poynting flux feeding 
the bubbles, and for astrophysical sources tb should be associated with a 
substantial variation in the outflow launching activity; observations of knots 
kinematics suggest characteristic times tb ~ 5 — 20 years; the experiments are 
in a similar regime . Because both time-scales and are relatively longer than 
the characteristic Keplerian period of rotation at the inner disk radius, jet 
launching should have ample time to reach steady-state. The characteristic 
astrophysical flow velocities can be taken to be w '--^ 200 km s~^. With these 
conditions, the presence of multiple bubble- like features should be observed on 
scales ranging from a few tens to a few hundred AU from the source. Indeed 
ejection variability, limb-brightened bubble-like structures and the presence of 
wiggles in the optical DC Tau jet are evident on scales ranging from of a few 
tens to a few hundred AU the source [H HH] . The experiments also indicate 
that asymmetries in the flow can be produced by instabilities that do not de- 
stroy the collimation, and because of their relatively short growth time , jets 
should develop non-axisymmetric features already within a few tens AU from 
the source, and become more heterogeneous and clumpy as they move further 
away to hundreds of AU. It was recently reported for a number of TTauri jets, 
including DC Tau, that already within 100 AU from the source the jet phys- 
ical conditions show considerable asymmetries with respect to the axis [14]. 
Finally over the same length scales the experiments suggest magnetic energy 
dissipation, heating of the plasma and a transition to a kinetically dominated 
jet which propagates ballistically. X-ray emission from the DC Tau jet was 
recently detected on the same length scales and it was proposed that mag- 
netic energy dissipation may be behind the heating mechanism [24] . As in the 
experiments, instabilities and the tangling of the magnetic field may provide 
a compelling route to efficient heating of such plasmas. 

4.3 Interaction with the interstellar medium 

Curved jets 

A number of bipolar Herbig-Haro (HH) jets exhibit a distinguishing C-shape 
morphology indicative of a steady bending [5] . Less regular curvature is also 
observed in a number of other HH jets; for example in HH 30 a small side 
drift close to the jet source is followed further away by a sudden bending 
[3]. In general, the curvature in jets has been linked either with the motion 
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of the jet sources relative to the ambient medium or with the presence of 
a widespread outflow; both cases giving rise to an effective transverse ram 
pressure (cross-wind) which curves the jet. Expected wind velocities vary from 
a few km for the jet-wind interaction associated with relative motions of 
TTauri stars with respect to the surrounding environment, (see for example 
|27j ) to typically higher velocities for irradiated jets, where best fits to HH505 
Ha emission maps were obtained for a wind velocity of 15 km s~^ |38j and 
estimates in [5] give wind velocities in the Orion nebula and in NGC1333 of 
-10-20 km s-i. 

The hydrodynamic laboratory jets described in section 4.1 are ideally 
suited to study the interaction with an ambient medium. The region into 
which the jet is launched (cf. Fig. [3| is a large vacuum that can be easily 
filled with different types of background gases. To investigate the dynamics 
of curved HH jets, a cross-wind was produced by a radiatively ablated foil 
appropriately placed in the jet propagation region [TJ 1321 130j . Typical wind 
velocities 30 — 50 km s^^ can be produced in the laboratory, with the impor- 
tant parameters characterising the interaction in the range Vjet/Vwind — 2 — 4 
and njf.t/n^ind ~ 0.1 — 10. Fig. 10 shows an example of the experiments and 



simulations of curved jets [5]. The characteristic dynamics of the interaction 
is similar for the laboratory and astrophysical systems, showing notably the 
formation of new working surface in the jet and "knotty" structure in the 
fiow. Curved jets are also Rayleigh- Taylor (RT) unstable, with the growth of 
such mode disrupting their propagation. However it was shown that jet ro- 
tation may partially suppress the instability by shearing the RT modes and 
confining the perturbations to a narrower layer of the jet body. Nevertheless 
this promotes the development of the Kelvin-Helmholtz instability (at least 
for the subsonically rotating jets) which is later responsible for disrupting 
the jet. Experimentally the RT growth time is of the order of the dynamical 
time over which the interaction can be produced and new, longer time-scale 
experiments will be needed to observe its full development. 



Clump propagation 

We now return to the scaling issue of some of the laboratory flows and in par- 
ticular the MHD jets. As we have seen their evolution is dominated by current- 
driven instabilities, and the resulting flow is inherently time-dependent and 
inhomogeneous. To study the propagation of such flows in an astrophysical 
setting, the data obtained from laboratory MHD jet simulations can be scaled- 
up and used as initial conditions to model astrophysical clumpy jets. There 
is clearly some arbitrariness on the choice of some of the scaling parameters 
and for the case presented here we assume the flow to be close to the YSO 
source. Noting that the laboratory and astrophysical jet velocities are of the 
same order, we choose the following three scaling: Viab = Vastro, Liab = 1 mm 
Lastro = 10 AU and put = 10"^ g cm'^ p^^t,,^ = lO'^* g cm'^. The 
choice of scale-length gives an initial jet radius ~20 - 30 AU and we take the 
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maximum jet density in the laboratory to scale to a maximum astrophysical 
jet density of ~10^ cm^^. The constraints on the scaling discussed in sec. [s] 
give: 1 ns ^ 0.05 years, 50 eV 3000 K, 50 T ^ 15 mG. The specific scahng 
applied in this case implies that the experimental flow, which lasts ~200 ns, 
corresponds to an astrophysical outflow lasting ~10 years. Such short times, 
when compared to the lifetimes of tens of thousands of years for protostel- 
lar outflows, may correspond to the ejection of a single "clump" as part of a 
more extended jet. The laboratory jet profiles to be scaled up are taken at 
a time approximately corresponding to the image in Fig. [TJd. The simulated 
astrophysical jet was evolved with the inclusion of cooling over ^50 years on 
a Cartesian grid of 400 x 10^ cells with a resolution of 2 AU. In these simu- 
lations the magnetic field is not included, and for the regime modelled here 
we would expect their inclusion to modify somewhat the dynamics. The flow 
dynamics is shown Fig. [TT] Initially the jet elongates because of the velocity 
variations imparted by the current-driven instability. The stretching of the 
jet is then followed by a rapid break up into smaller clumps which move at 
different velocities. The structure appearing in the large knot forming in the 
flow appears to be the result of Rayleigh- Taylor instabilities. In general, a sin- 
gle dense clump is produced by an ejection event like this, with the resulting 
outflow remaining well coUimated over the propagation across ^3000 AU. 

5 Summary 

Progress in high-energy density plasma experiments on lasers and z-pinch fa- 
cilities has permitted in the last ten years to start investigating a range of 
"large-scale" astrophysical phenomena in the laboratory; extending the tradi- 
tional domain of laboratory astrophysics beyond the work on micro-physics. 
Through careful design of the experiments, the plasma produced can be scaled 
to the astrophysical environment; allowing complex, intrinsically non- linear, 
three-dimensional phenomena to be accessed in a controlled manner. An im- 
portant outcome being the validation of astrophysical codes on the labora- 
tory data. Although work on astrophysical jets has been performed on both 
lasers and z-pinch facilities, we have focused here only on the studies of jets 
produced on the MAGPIE z-pinch facility. Two main "types" of jets were 
developed: hydrodynamic jets, to be used for propagation studies, and mag- 
netohydrodynamic jets of interest to the launching phase. In general, there 
is some considerable control on the experiments: the initial condition can be 
partially modified, for example the density and magnetic field distributions; 
more complex physics, such as rotation, can be added, and different plasma 
condition can be produced, by modifying for example the cooling rates. Over- 
all the combination of laboratory experiments and simulations can provide 
some very important insights on the physics of astrophysics, and as technol- 
ogy advances we can expect evermore exotic phenomena to be reproduced in 
the laboratory. 
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Table 1. Characteristic conditions in laboratory (Z-pinch) and YSO jets 
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Fig. 1. Plot of the temperature versus density for a variety of laboratory and astro- 
physical plasmas. Lines of constant pressure are shown for fully ionized hydrogen. 
High-energy density laboratory astrophysics is in the regime of pressures > 1 Mbar 
and the typical conditions currently obtained in laser and z-pinch experiments are 
easily in this range. The conditions that will be accessible in future Magnetic Fusion 
(MF) reactors and Inertial Confinement Fusion (ICF) laser experiments are also 
indicated. The main phases of the Interstellar Medium (ISM) are shown. The Inter- 
galactic Medium (IGM) and the Intracluster Medium (ICM) lie outside the plot at 
lower densities. 
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Fig. 2. The lower part of each plot shows schematically the "engine" that pro- 
duces astrophysical (left) and laboratory jets (right). The physical mechanisms are 
clearly very different, however we are interested in the scaling of the flows inside the 
Modelling Boxes. 
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Fig. 3. (a) Typical arrays are made with 16 tungsten wires, each with a diameter of 
18 ^m. The smaller radius of the array is 8 mm and the wires are inclined at an angle 
of 30° with respect to the axis. The axial length of the array is 12 mm. Continuous 
plasma streams converge on axis producing a conical shock which redirects the flow 
axially. (b) Time-resolved, filtered XUV emission from a laboratory jet. Emission 
from the region in between the dotted lines is screened by the anode. As the jet 
propagates it cools down and its emission decays, however it is again visible where 
the jet interacts with a background plasma. 
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Fig. 4. Tiie experimental set-up is similar to that shown in Fig. [3| in this case 
however the wires are twisted in the azimuthal direction (see main text) (a) Az- 
imuthally averaged profiles of azimuthal velocity, density and electron temperature 
as a function of radial position, taken 6 mm above the cathode 200 ns after the 
start of the current pulse. The profiles are for different twist angles of the array, 
which can be used to modify the angular rotation in the jet. (b) Isodensity surfaces 
at 270 ns show the dense plasma around the wires (10~^ g cm , dark gray) and 
plasma streams (2.5 x 10"^ g cm , light gray); velocity streamlines from one of the 
wires are shown with the oppositely rotating flows separated visually by red and 
orange streamlines. The oppositely rotating flow is due to electrode's effects and it 
is dynamically unimportant. The azimuthally averaged poloidal magnetic field lines 
are shown in blue.|2] 
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Fig. 5. Schematic of a radial wire array. 
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Fig. 6. Slices of mass density from a 3D MHD simulation of radial wire arrays. The 
evolution is shown at times (a) 165 ns, (b) 175 ns, (c) 185 ns and (d) 205 ns. The 

logarithmic density scale is from 10~^ g cm (blue) to 10~^ g cm (red). Regions 
in white are void of plasma, but not electromagnetic fields, and are essentially a 
computational "vacuum". The square, black regions are the electrodes. |11| 
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Fig. 7. (a) Magnetic field lines (yellow) and current density (red) distribution inside 
the magnetic cavity at 245 ns To show the inside of the magnetic cavity the 
isodensity contours (shades of blue) are sliced vertically. With the onset of the 
kink instability the magnetic field wraps tightly around the jet, which is seen more 
clearly in (b) without the cavity walls [12j. (c) Experimental shadowgraph showing 
the bubble wall and the clumpy jet launched along the axis|33|. 
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Fig. 8. 3D MHD simulation of a radial foil experiment. The images show the line- 
of-sight integrated emission (in arbitrary units) at different times. The red line in 
the first panel shows the position of the electrodes. 




Fig. 9. Experimental image sliowing the presence of nested magnetic cavities with 
an embedded jet. The image shows the self-emission in the XUV range [13] 
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Fig. 10. (a) Experimental time-resolved XUV image of a curved jet. (b) Synthetic 
XUV image of a curved laboratory jet from a 3D simulation, (c) Column density 
from a 3D simulation of an astrophysical jet[S] 
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Fig. 11. Propagation of an episodic protostellar jet showing the break up into small 
and large clumps. 



